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Dynamic fluctuations in the local density of non-identified hadron tracks re- 
constructed in the STAR TPC are studied using the discrete wavelet transform 
power spectrum technique which involves mixed event reference sample compari- 
son. The two-dimensional event-by-event analysis is performed in pseudo-rapidity 
rj and azimuthal angle <fi. HIJING simulations indicate that jets and mini-jets 
result in signals, visible without high px selection, when the dynamic texture 
analysis is applied. Scanning a broad range of event multiplicities, we study the 
dependence of the signals on the initial conditions. Event structures are studied 
separately with positive and negative tracks, as well as both charges. A change of 
regime is observed in AuAu collisions at v5aw = 130 GeV as event multiplicity 
is increased: a long range rj correlation (or suppressed fluctuation vis-a-vis mixed 
events) is seen in same charge data. This effect is qualitatively similar to one of 
the predicted manifestations of the Color Glass Condensate. 

1 Introduction 

Bulk properties of strongly interacting matter under extreme conditions are 
the focus of the on-going RHIC program. Deconfmement and chiral symmetry 
restoration []J are expected to take place in collisions of ultra- relativistic nuclei. 
Because these phase transitions are multiparticle phenomena, a promising, 
albeit challenging, approach is the study of dynamics of large groups of final 
state particles. The dynamics shows itself in the correlations and fluctuations 
(texture) on a variety of distance scales in momentum space. 

The multi-resolution dynamic texture approach (applied for the first time|| 
at SPS) uses discrete wavelet transform P](DWT) to extract such information. 
At present stage, the information is extracted in a comprehensive way, without 
any built-in assumptions or filters. Mixed events are used as a reference for 
comparison in search for dynamic effects. Event generators are used to "train 
intuition" in recognizing manifestations of familiar physics (such as elliptic 
flow or jets) in the analysis output, as well as to quantify sensitivity to the 
effects yet unidentified, such as critical fluctuations or clustering of new phase 
at hadronization. 




Figure 1: Haar wavelet basis in two dimensions. The three modes of direc- 
tional sensitivity are: a) diagonal b) azimuthal c) pseudo-rapidity. For the 
finest scale used, the white rectangle drawn "on top" of the function in panel 
a) would correspond to the smallest acceptance bin (pixel). Every subsequent 
coarser scale is obtained by expanding the functions of the previous scale by 
a factor of 2 in both dimensions. (Reproduced from |2|). 

2 The STAR experiment 

STAR Time Projection Chamber |4]] (TPC), mounted inside a solenoidal mag- 
net, tracks charged particles within a large acceptance (\r)\ < 1.3, < < 2ti) 
and is well suited for event-by-event physics and in-depth studies of event 
structure. The data being reported are obtained during the first (\/Snn = 130 
GeV) year of RHIC operation. The minimum bias trigger discriminates on a 
neutral spectator signal in the Zero Degree Calorimeters 0. By adding a re- 
quirement of high charged multiplicity within |?7| < 1 from the scintillating 
Central Trigger Barrel, one obtains the central trigger. Vertex reconstruction 
is based on the TPC tracking. Only high quality tracks found to pass within 

3 cm of the event vertex are accepted for the texture analysis. 

3 Dynamic texture analysis procedure 

Discrete wavelets are a set of functions, each having a proper width, or scale, 
and a proper location so that the function differs from only within that 
width and around that location. The set of possible scales and locations is 
discrete. The DWT transforms the collision event in pseudo-rapidity 77 and 
azimuthal angle <fi into a set of two-dimensional functions. The basis functions 
are defined in the (77, 0) space and are orthogonal with respect to scale and 
location. We accumulate texture information by averaging the power spectra 
of many events. 



The simplest DWT basis is the Haar wavelet, built upon the scaling func- 
tion g{x) = 1 for < x < 1 and otherwise. The function 

f{x) = {+1 for < x < — ; — 1 for — < x < 1; otherwise} (1) 

is the wavelet function. 

The experimental acceptance in 77, 0, and px ((M < 1, < < 2n)) is 
partitioned into bins. The r/-0 partitions are of equal size, whereas in p?, the 
binning is exponential when more than one pt bin is used. In each bin, the 
number of reconstructed tracks satisfying the quality cuts is counted. 

The scaling function of the Haar basis in two dimensions (2D) G(<f), rj) = 
g(4>)g(r]) is just a bin's acceptance (modulo units). The wavelet functions 
F x (where the mode of directional sensitivity A can be azimuthal 0, pseudo- 
rapidity r/, or diagonal (prj) are 

F* n = mm, ** = f&Mv), F* = gWiv). (2) 

We set up a two dimensional (2D) wavelet basis: 

< M .(0, V ) = 2 m F A (2"> - i, 2 m V - j), (3) 

where m is the integer scale fineness index, and % and j index the positions of 
bin centers in and r/. Then, i • with integer m, i, and j are known || 
to form a complete orthonormal basis in the space of all measurable functions 
defined on the continuum of real numbers L 2 (M). We construct G miij (0, 77) 
analogously to Eq.|3|. 

Fig. [I] shows the wavelet basis functions F in two dimensions. At first 
glance it might seem surprising that, unlike the ID case, both / and g enter the 
wavelet basis in 2D. Fig. [I] clarifies this: in order to fully encode an arbitrary 
shape of a measurable 2D function, one considers it as an addition of a change 
along (f(4>)g(ri), panel (b)), a change along rj (c/ (0) / (77) , panel (c)), and 
a saddle-point pattern {f{4>)f{f])i panel (a)), added with appropriate weight 
(positive, negative or zero), for a variety of scales. The finest scale available is 
limited by the two track resolution, and, due to the needs of event mixing, by 
the number of available events. The coarser scales correspond to successively 
re-binning the track distribution. The analysis is best visualized by considering 
the scaling function G TOi ij(0, rj) as binning the track distribution p(0, rf) in bins 
i,j of fineness m, while the set of wavelet functions F^ li j(<j), rf) (or, to be exact, 
the wavelet expansion coefficients (p, F^^)) gives the difference distribution 



between the data binned with given coarseness and that with binning one step 
finer. We use WAILI|| software to obtain the wavelet expansions. 

In two dimensions, it is informative to present the three modes of a power 
spectrum with different directions of sensitivity P^ v (m), P^(m), P v (m) sepa- 
rately. We define the power spectrum as 

where the denominator gives the meaning of spectral density to the observ- 
able. So defined, the P x (m) of a random white noise field is independent of 
m. However, for physical events one finds P x {m) to be dependent on m due 
to the presence of static texture features such as acceptance asymmetries 
and imperfections (albeit minor in STAR), and non-uniformity of the dNj dr\ 
shape. In order to extract the dynamic signal, we use P x (m) true — P x (m) m i X 
where the latter denotes power spectrum obtained from the mixed events. 
The mixed events are composed of the (77, <fi) pixels of true events, so that a 
pixel is an acceptance element of the finest scale used in the analysis, and in 
no mixed event is there more than one pixel from any given true event. The 
minimum granularity used in the analysis is 16 x 16 pixels. Q 

Systematic errors can be induced on P x (m) true — P x (m) mix by the process 
of event mixing. For example, in events with different vertex position along 
the beam axis, same values of r\ may correspond to different parts of the TPC 
with different tracking efficiency That will fake a dynamic texture effect in 77. 
In order to minimize such errors, events are classified into event classes with 
similar multiplicity and vertex position. Event mixing is done and P x (m) true — 
P x ijn)mix is constructed within such classes. Only events with z vertex lying 
on the beam axis within 30 cm from the center of the chamber are accepted for 
analysis. To form event classes, this interval is further subdivided into five bins. 
We also avoid mixing of events with largely different multiplicity. Therefore, 
another dimension of the event class definition is that of the multiplicity of 
quality tracks in the TPC. For central trigger events, the multiplicity range 
of an event class is typically 25. For small multiplicity events taken with 
the minimum bias trigger, the multiplicity range per event class is 10 for 
multiplicities below 90, and 20 in the range between 90 and 150. For larger 
multiplicities in the minimum bias sample, this multiplicity range was taken to 
be 50. Events with less than 15 good quality tracks in the fiducial acceptance 
are ignored. 

1 For a quick reference, here are the scales in -q. Scale 1: Arj — 1; scale 2: A77 = 1/2; 
scale 3: A77 = 1/4 and so on. 
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Figure 2: (P^ rue — Pmix) fornix (scale 1) in Au Au collisions at y/s = 130 GeV 
as a function of normalized multiplicity for all charged particles. * - STAR; 
• - regular HIJING; o - HIJING without jets. 

4 Dynamic textures in the STAR data 

Elliptic flow is a prominent large scale dynamic texture effect already well 
measured at RHIC0. The DWT approach localizes elliptic flow on scales 2 
and, to some degree, 3 of the Pf rue — P^ ix - In this report, we ignore flow and 
concentrate on the r\ observables. 

Multiplicity scans reveal dependence of the signals on the initial conditions: 
number of participants and binary collisions and the energy density and size 
of the interacting system. Such scans as a function of normalized multiplicity 
are shown in Fig. |2|, ^, and |J The figures show the scale 1 (large scale) 
observables for tracks of both charges, and for the positive only and the nega- 
tive only, respectively. The horizontal error bars on the points are formed by 
multiplicity boundaries of event classes used in the analysis or, when a coarser 
re-binning of the multiplicity is done on the analyzed data to enhance the 
presentation, reflect that re-binning. The event mixing is always confined to 

2 To form the normalized multiplicity, the multiplicity n of quality tracks is divided by 
??o, where no is such that 99% of minimum bias events have multiplicity of quality tracks 
less than no . The normalized multiplicity depends weakly on the experimental definition of 
minimum bias and on the track quality cuts. 




Figure 3: (Pf rue — Pmix)/Pmix (scale 1) in Au+Au collisions at s/s = 130 GeV 
as a function of normalized multiplicity for positively charged particles. * - 
STAR; • - regular HIJING; o - HIJING without jets. 



the specific multiplicity and z-vertex position class, described in Section y. 

A comparison Q with HIJING|| is done in order to understand the effects 
of energy-momentum conservation, jets and resonances on our measurements. 
For the regular HIJING and HIJING without jets, the multiplicity dependence 
is weak or absent. This is understandable given the nature of HIJING as a 
super-position of binary nucleon-nucleon collisions with no collective effects. 
The regular HIJING and HIJING without jets show different magnitude of 
the signal. This underscores the role of jets in creating local density fluctua- 
tions (positive dynamic texture) in HIJING. Notice that this result is obtained 
without imposing a high cut. The experimental data are close to regular 
HIJING for peripheral events, but their texture becomes suppressed as multi- 
plicity grows. 

To understand the nature of this suppression, we look at same charge 
particles (Fig. ^|and^). We see that the STAR data are not really intermediate 
between the regular HIJING and HIJING without jets, but display a more 

3 In the event generator data in the figures, no GEANT and no response simulation 
is done. Instead, only stable charged particles (e,fj,,ir,K ,p) and their antiparticlcs from 
the generator output are considered, provided that they fit into the STAR TPC fiducial rj 
acceptance |ry| < 1. Momentum resolution and pt acceptance arc not simulated. 



: STAR prelir 

L T 1 I T T t 
-Til **^u 


ninary 

it 

T °T 


O 1 i— I— i 
" 


#^##IT — 

^ t 1 I tWitP t 1 -U 
I 4JLt TV i 1 

i— A— i ^ 
1 . . . 1 



0.2 0.4 0.6 0.8 1 

multiplicity/!!,, 



Figure 4: (P^ ue — Pmix)/Pmix (scale 1) in Au+Au collisions at yfs = !30GeV 
as a function of normalized multiplicity for negatively charged particles. * - 
STAR; • - regular HIJING; o - HIJING without jets. 



complex behavior. This is seen from the fact that as the multiplicity grows, the 
dynamic texture becomes negative and stays roughly constant past a transition 
region which occupies about a half of the multiplicity range, whereas the "no 
jets" prediction points are consistent with for all multiplicities. How to 
explain the negative dynamic texture? Recall the definition of the observable 
(Section |3|). At scale 1, we are looking at the fluctuation in same-charge 
track occupancy between two "pixels", each one unit of rj long, located in 
forward/backward hemispheres around mid-rapidity. The negative Pt rue —PZix 
means that in the mixed events, this occupancy fluctuation is stronger, i.e. a 
correlation takes place in the real events. 

4 This looked very counter-intuitive to us at first - indeed, the events are composed out of 
positive and negative charge sub-events. If the positive sub-event shows reduced (negative) 
fluctuations, and so does the negative sub-event, how can the fluctuation measure for the 
combined event be above zero ? To ensure that the "balance" is kept, we construct and 
measure from p and p m i X for positive and negative tracks. In the following, the (p\F) 
notation denotes coefficients in the expansion of an individual event track density p into 
wavelet basis F. For simplicity, we omit the indices. We sum over locations but not over 
the scale index. Averaging over many events is assumed everywhere. 

(p+ + p-\F) 2 - (p+ lx + Pmlx \F) 2 = 



The tracking conditions in the chamber are different at high and low mul- 
tiplicities, and the effect of that on the measurements in question needs to be 
understood. We do this by processing HIJING events through full GEANT 
and response simulation, and applying the actual experimental reconstruction 
to those events. Such processing (done for a smaller number of high multiplic- 
ity HIJING events than is shown in the figures) reveals no significant difference 
for the reported observable. 

Thus, we are looking at an interplay of fluctuation and correlation effects, 
neither of which is trivial, as a function of initial conditions. The signals 
become weaker for finer scales. At this point one can only speculate about the 
nature of the correlation effect in central events, but the HIJING simulation 
indicates that jet quenching (or even total disappearance of jets) does not 
account for the negative Pj rMe — P m i X , even though it may be a prerequisite 
for its observation, given that jets work to create a positive texture. What 
about Bose-Einstein correlation? In the simulations, its effect can not be 
introduced into the multiparticle distributions event-by-event from the first 
principles, and in HIJING it is not considered at all. Low values of Q% nv (< 50 
MeV@), typical for both Bose-Einstein and Coulomb effects, make it unlikely 
for these effects to be responsible for correlations with characteristic A77 « 1. 

Longitudinal expansion maps rapidity differences onto time differences; the 
features characterized by larger rapidity differences must be formed early. 
In the Color Glass Condensate picture, multiplicity fluctuations in differ- 
ent rapidity windows are predicted to be correlated for large (I /as) rapidity 



intervals [10]. This is seen as a consequence of the classical coherence of the 
gluon field. Such or a qualitatively similar effect is indeed required to explain 
the data. 



{p + \F) 2 + (p-\Ff + 2{p+\F){p-\F) (p+ llx \Ff - (p- ix \Ff - 2{p+ llx \F) { Pmlx \F) (5) 



We measure independently (p + \F) — (p^, lix \F) (Fig|| modulo normalization) and 
(p~\F) 2 — (p mix \F) 2 (Fig.|], modulo normalization). In order to test that EqJ| holds, 
we need to obtain the correlation terms between the opposite charges (p + \F)(p~\F) and 
(Pmix\F) (Pmix\F) ■ This can be done by comparing power spectra of images filled with equal 
weights for positive and negative particles, and of those where negative particles are entered 
with a negative weight: 

(p+ + p-\F) 2 (p+ p-\F) 2 = 4(p+\F)(p-\F), (6) 

which is true for p m ix as well. With this input, validity of EqJ| has been established and 
the "paradox" presented above resolved: it is the large scale correlation between positive 
and negative charges that accounts for the "extra" texture when the positive and negative 
sub-events arc combined. 



5 Conclusions 

The STAR measurements of AuAu data reveal a non-trivial picture of non- 
statistical correlations and fluctuations which is qualitatively different for pe- 
ripheral and central collisions. A possible interpretation of the central data is 
suppression of jets, combined with particle emission correlated over the length 
of the order of a unit of rapidity. 
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